Narrow Line Photoassociation in an Optical Lattice 
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With ultracold 88 Sr in a ID magic wavelength optical lattice, we performed narrow line photoas- 
sociation spectroscopy near the 1 So— 3 Pi intercombination transition. Nine least-bound vibrational 
molecular levels associated with the long-range 0„ and l u potential energy surfaces were measured 
and identified. A simple theoretical model accurately describes the level positions and treats the 
effects of the lattice confinement on the line shapes. The measured resonance strengths show that 
optical tuning of the ground state scattering length should be possible without significant atom loss. 

PACS numbers: 34.80.Qb, 32.80-t, 32.80.Cy, 32.80.Pj 



Photoassociation (PA) spectroscopy [jj, |2( is a valu- 
able tool for studies of atomic collisions and long-range 
molecules 0, 0, E3, 0, II 0. In contrast to con- 
ventional molecular spectroscopy that probes the most 
deeply bound vibrational levels of different electronic 
states, PA spectroscopy excites the vibrational levels 
close to the dissociation limit, where the molecular prop- 
erties depend mostly on the long-range dipole-dipole in- 
teratomic interactions. Molecular potentials are much 
simpler at long range than at short range, and are di- 
rectly related to the basic atomic properties such as the 
excited state lifetime and the s-wave scattering length. 

In this Letter, we report the use of PA spectroscopy 
to resolve nine least-bound vibrational levels of the 88 Sr2 
dimer near the 1 So— 3 Pi intercombination line. Our iden- 
tification of these levels is required for tuning of the 
ground state scattering length with low-loss optical Fes- 
hbach resonances [ljj. This is of great interest for Sr, 
since magnetic Feshbach resonances are absent for the 
ground state, and the background scattering length is 
too small to allow evaporative cooling || . In contrast to 
prior PA work that utilizes strongly allowed transitions 
with typical line widths in the MHz range, here the spin- 
forbidden 1 So— 3 Pi line has a natural width of 7.5 kHz. 
This narrow width allows us to measure the least-bound 
vibrational levels that would otherwise be obscured by a 
broad atomic line, and to observe characteristic thermal 
line shapes even at /iK atom temperatures. It also per- 
mits examination of the unique crossover regime between 
the van der Waals and dipole-dipole interactions, that 
occurs in Sr near the 1 So— 3 Pi dissociation limit. This 
access to the van der Waals interactions ensures large 
bound-bound Franck-Condon factors, and may lead to 
more efficient creation of cold ground state Sr2 molecules 
with two-color PA than what is possible using broad tran- 
sitions. Bosonic 88 Sr is ideal for narrow line PA because 
of its hypcrfmc-frcc structure and a high isotopic abun- 
dance (83%). In addition, Sr is a candidate atom for op- 
tical clocks [ill 0] , and PA spectroscopy could help 
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FIG. 1: Schematic diagram of the long-range Sr2 molecular 
potentials (not to scale). The ground and excited molecu- 
lar states coincide with two ground state atoms, and with 
one ground and one excited state atom, respectively, at large 
internuclear separations. The ground state has gerade sym- 
metry and its energy is given by the potential V g , while the 
excited state ungerade potentials that support dipole tran- 
sitions to the ground state are Vou and Vi u , the latter with 
a small repulsive barrier. All vibrational states of U and 
lu (dashed and dotted lines, respectively) are separated by 
more than the natural line width, permitting high resolution 
PA spectroscopy very close to the dissociation limit when the 
atoms are sufficiently cold. 



determine density shifts for various clock transitions. 

To enable long interrogation times needed for narrow 
line photoassociation, we trap Sr atoms in a ID opti- 
cal lattice. This technique also suppresses recoil shifts 
and Doppler broadening and thus simplifies the PA line 
shapes. The specially chosen lattice wavelength (A mag i c ) 
induces state-insensitive AC Stark shifts that minimize 
the perturbing effects of the lattice Q . 

FigureQ]illustrates the relevant potential energy curves 
for the Sr2 dimer as a function of interatomic separation 
R. The photoassociation laser at 689 nm induces allowed 
transitions between the separated 1 So atom continuum at 
a temperature T to the vibrational bound levels of the 
excited ungerade potentials Vq u and V\ u , corresponding 
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FIG. 2: Experimental configuration. The 2 /iK atoms are 
confined in a ID optical lattice in the Lamb-Dicke regime 
that allows Doppler- and recoil-free photoassociation spec- 
troscopy. The lattice laser intensity is large enough to confine 
most atoms in the ground motional state. The 689 nm weak 
PA laser co-propagates with the lattice. The 914 nm lattice 
wavelength ensures that So and the mj = ±1 sublevels of 
3 Pi are Stark-shifted by equal amounts. 

to the total atomic angular momentum projections onto 
the internuclear axis of and 1, respectively. Although 
our calculations, based on the coupled channels method 
in Ref. 0], account for the Coriolis mixing between the 
U and l u states, this mixing is small enough that U and 
l u are good symmetry labels for classifying the states. 
Including the Coriolis mixing is necessary in order to fit 
the observed 0„ and l u bound state energies. Only s- 
waves (total ground state molecular angular momentum 
J = 0) contribute to the absorption (excited state J = 1) 
at our low collision energies. The long-range potentials 
are given by the Cq /R 6 van der Waals and C3 / R 3 dipolc- 
dipole interactions, plus a rotational term cx 1/R 2 , 

V 0u =-C 6 , 0u /R 6 -2C 3 /R 3 + h 2 A 0u /(8Tr 2 iiR 2 ), (1) 
Viu =-C\ lu /R 6 + C 3 /R 3 + h 2 A lu /(8n 2 fiR 2 ), (2) 

where h is the Planck constant, /i is the reduced mass, 
Aq u = J(J + 1) + 2, and A lu = J(J + 1). The values of 
C3, Ceflu, and Cq^ u are adjusted in our theoretical model 
so that bound states exist at the experimentally deter- 
mined resonance energies. Since our modeling of the PA 
energy levels due to long-range interactions is not sensi- 
tive to the details of short-range interactions, the latter 
are represented by simple 6-12 Lennard- Jones potentials 
whose depths are chosen to approximately match those 
in Ref. |lj| ■ The coefficient C3 can be expressed in terms 
of the atomic lifetime r as C3 = 3fic 3 / (Atlu 3 ) , where Tllo 
is the atomic transition energy and c is the speed of light. 
Since the C3 and C§ interactions have comparable mag- 
nitudes in the region of R relevant to the energy levels 
we observe, the standard semiclassical analysis based on 
a single interaction cx R~ n 17] is not applicable. 

The experiment is performed with ~ 2 /iK 88 Sr atoms 
that are trapped and cooled in a two-stage magneto- 
optical trap (MOT) [Uli^. As the atoms are cooled in 
the 689 nm narrow line 1 So— 3 Pi MOT, they are loaded 
into a ID far-detuned optical lattice |l3(, which is a 300 
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FIG. 3: The spectra of the long-range Sr2 molecule near the 
1 So— 3 Pi dissociation limit. The horizontal scale is marked 
on the rightmost panel; different PA laser intensities were 
used for each spectrum due to largely varying line strengths. 
The top labels indicate the interatomic separations that corre- 
spond to the classical outer turning points of each resonance. 



mW standing wave of 914 nm light with a 70 /im beam 
waist at the MOT. This results in 10 5 atoms with the 
average density of about 3 x 10 12 /cm 3 , after the MOT 
is switched off. The axial trapping frequency in the lat- 
tice is v z ~ 50 kHz, much larger than the 5 kHz atom 
recoil frequency and the 7.5 kHz line width, resulting in 
Doppler- and recoil-free spectroscopy if the PA laser is 
collinear with the lattice beam. Figure |2] illustrates the 
lattice configuration. The polarizations of the lattice and 
PA lasers are linear and orthogonal in order to satisfy the 
Amagic condition at 914 nm |l4j. The choice of A mag i c en- 
sures that 1 So and the magnetic sublevels of 3 Pi that are 
coupled to 1 So by the PA laser are Stark shifted by equal 
amounts, resulting in the zero net Stark shift for the tran- 
sition and consequently eliminating inhomogeneous Stark 
broadening and shifts of the PA lines. A 689 nm diode 
laser used for PA is offset-locked to a cavity-stabilized 
master laser with a sub-kHz spectral width [111 Il3| . 

To trace out the molecular spectra, the PA laser fre- 
quency is stepped, and after 320 ms of photoassocia- 
tion at a fixed frequency the atoms are released from 
the lattice and illuminated with a 461 nm light pulse 
(resonant with 1 So— 1 Pi) for atom counting. At a PA 
resonance, the atom number drops as excited molecules 
form and subsequently decay to ground state molecules 
in high vibrational states or hot atoms that cannot re- 
main trapped. Figure shows the nine observed PA line 
spectra near the dissociation limit. 

The line fits for the free-bound transitions are based on 
a convolution of a Lorentzian profile with a thermal dis- 
tribution of initial kinetic energies [2(|, as schematically 
shown in Fig. 01(a). In the ID optical lattice, we work in 
the quasi-2D scattering regime, since T ~ 2 /iK ~ v z /kB 
(ks is the Boltzmann constant), and the axial motion is 
quantized. For weak PA laser intensities used here, the 
total trap loss rate is a superposition of loss rates K e for 
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atoms with thermal energies e = h 2 k 2 / (8ir 2 fj) , 

Ke{v) = Ji(e/h + v-l) 2 + ( 7 + 7i) 2 /4' (3) 

where v and v§ are the laser and molecular resonance 
frequencies expressed as detunings from the atomic 
1 So— 3 Pi line, 7 ~ 15 kHz is the line width of the excited 
molecular state, 71 ~ 25 kHz is a phenomenological pa- 
rameter that accounts for the observed line broadening, 
l op t is the optical length 01 discussed below, and only 
s-wave collisions are assumed to take place. While the 
total 3D PA loss rate is 

the 2D loss rate is 

f°° rip 

Jo k B T 

Dimensional effects appear in Eq. © as a larger density 
of states at small thermal energies {i.e. no y/e factor), and 
as a red shift by the zero-point confinement frequency 
Ug/2. Another property of a 2D system (not observed 
here due to small probe intensities) is a lower bound on 
power broadening near T ~ 0, fixed by the zero-point 
momentum |2 1| . Equation JSJ) assumes only collisions in 
the motional ground state of the lattice potential. Al- 
though 50% of our atoms are in excited motional states, 
using a more complete expression yielded essentially the 
same line shapes. If any local atom density variations 
associated with PA are neglected and the temperature is 
assumed to be constant during probing, then the atom 
density evolution is given by 

h = -2Kn 2 - n/n, (6) 

where 77 ~ 1 s is the lattice lifetime, from which we 
calculate the number of atoms remaining in the lattice. 

The measured PA resonances Vq were reproduced the- 
oretically by adjusting the C3 and C$ parameters as well 
as the short-range interatomic potentials. The measured 
and calculated vq are compared in TableQ] The positions 



On Measured 


U Calculated 


1„ Measured 


l u Calculated 


(MHz) 


(MHz) 


(MHz) 


(MHz) 


-0.435(37) 


-0.418 


-353.236(35) 


-353.232 


-23.932(33) 


-23.927 


-2683.722(32) 


-2683.729 


-222.161(35) 


-222.152 


-8200.163(39) 


-8112.987 


-1084.093(33) 


-1084.091 






-3463.280(33) 


-3463.296 






-8429.650(42) 


-8420.124 







TABLE I: The comparison of measured and calculated molec- 
ular resonance frequencies, for both U and l u potentials. 



of the measured lines are found by fitting Eqs. (l.'iKit to 
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FIG. 4: (a) A line is schematically shown as a sum of 
Lorentzians with positions and amplitudes given by a con- 
tinuous thermal distribution, (b) The measured line shapes 
of the —222 MHz transition at 1.7 /iK (open circles) and 3.5 
fiK (filled circles) , clearly showing the effect of temperature on 
the red side of the line. The bottom trace shows the residuals 
of the 1.7 /iK curve fit. 



the data (see Fig. |2J|. Both 2D and 3D line shape mod- 
els were applied, without an appreciable difference in the 
obtained optical lengths. However, the 2D model yields 
Vq corrections of +15(10) kHz that are included in Table 
[I] The measurement errors in the Table are a combi- 
nation of a 30 kHz uncertainty of the atomic reference 
for the PA laser, and of the systematic and statistical 
reproducibility. The calculated vq results are based on 
the coupled channels model 15] and were obtained us- 
ing C 3 = 0.007576 au (r = 21.35 /is), C 6 . 0u = 3550 au, 
and Cq.i u — 3814 au (1 au is E^a^ for C3 and E^a® 
for C 6 , where E h = 4.36 x 10~ 18 J and a = 0.0529 
nm). However, the model limitations, such as insuf- 
ficient knowledge of the short-range potentials, require 
the following uncertainties in the long-range parameters: 
C 3 = 0.0076(1) au (r = 21.4(2) /its), C 6:0u = 3600(200) 
au, and Ce,i u — 3800(200) au. The theoretical results for 
the two PA lines at detunings over 8 GHz disagree with 
experiment by 0.1-1% due to a high sensitivity to the 
short-range potentials, while all the other lines fit well 
within the specified uncertainties of as low as 10~ 5 . 

Since the width of the Sr intercombination line is in 
the kHz range, even ultracold temperatures of a few /iK 
contribute significant thermal broadening. Figure ^ (b) 
shows two spectra of the —222 MHz PA line taken at 
T = 1.7 /tK and T — 3.5 /iK (T was measured from the 
time-of- flight atom cloud expansion). The spectrum of 
the hotter sample manifests thermal broadening as a tail 
on the red-detuned side. 

A photoassociation line depth is given by the optical 
length 1Q], Z op t, a parameter proportional to the free- 
bound Franck-Condon factor and the PA laser intensity. 
It is useful for specifying the laser-induced changes in 
both elastic and inelastic collision rates that lead to op- 
tical control of the atomic scattering length and to atom 
loss, respectively. The optical lengths were obtained from 
the data and calculated from the Franck-Condon factors 
under the assumption of a zero ground state background 
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scattering length |9j. For the —0.4 MHz line, the ex- 
perimental lopt/I — 4.5 x 10 5 flo/(W/cm 2 ), while the 
theoretical value is 9.0 x 10 5 ao/(W/cm 2 ), where / is the 
PA laser intensity. Significant experimental uncertainties 
may arise from the calibration of the atomic density and 
the focused PA laser intensity. In addition, we observe 
line broadening 71, most likely a result of the residual 
magnetic field, and any variation in this extra width can 
lead to Z op t error. The —0.4 MHz resonance is special 
due to its very large classical turning point, and its ef- 
fective optical length increases with decreasing e, while 



opt is energy-independent for the other lines (due to the 



Wigner threshold law |22j). Therefore, the reported ef- 
fective Z op t values for the —0.4 MHz line are based on 
the 2 [iK sample temperature with the 1 /xK zero-point 
energy. In addition, the fitted Z op t for this resonance was 
multiplied by 2.5 to obtain the value quoted above, since 
our calculations show that 60% of the photoassociated 
and decayed atoms have insufficient kinetic energies to 
leave the ~ 10 /iK deep trap. The measured and cal- 
culated Z pt for the other PA lines can differ by up to 
an order of magnitude, but both decrease rapidly as the 
magnitude of the detuning increases and are of the order 
of Z op t /I = 1 ao/(W/cm 2 ) for the most deeply bound lev- 
els. The theoretical optical length of the —222 MHz line 
was found to have the largest dependence on the assumed 
ground state scattering length due to its sensitivity to the 
nodal structure of the ground state wave function. 

Intercombination transitions of alkaline earths such as 
Sr are particularly good candidates for optical control of 
the ground state scattering length, a opt , because there is 
a possibility of large gains in a op t with small atom losses. 
In fact, using the —0.4 MHz PA line with the measured 
Z op t ~ 5 x 10 5 ao cm 2 /W will allow tuning the ground 
state scattering length by 0] a opt = Z op t(7/<5)/ — ±300 
ao, where the PA laser with / = 10 W/cm 2 is far-detuned 
by 8 — ±160 MHz from the molecular resonance, and 
the factor / = (1 + {i/25) 2 {\ + 2fcZ opt ) 2 )- 1 ~ 0.8 ac- 
counts for the power broadening for 3 /iK collisions. In 
contrast, optical tuning of the scattering length in alkali 
87 Rb |23j achieved a op t = ±90 ao at much larger PA laser 
intensities of 500 W/cm 2 . In addition, the Sr system at 
the given parameter values will have a loss rate of 01 
K = (2/i//i)a op t7/(2<J) ~ 2 x 10~ 14 cm 3 /s, while the loss 
rate in the 87 Rb experiment was 2 x 10~ 10 cm 3 /s. The 
overall optical length gain of over 5 orders of magnitude 
is possible for the Sr system because the narrow inter- 
combination transition allows access to the least-bound 
molecular state, and the PA line strength increases with 
decreasin g de tuning from the atomic resonance (see Fig. 
1 in Ref. [ill). The above a op t and K values accessible 
with the —0.4 MHz resonance result in the clastic and in- 
elastic collision rates of r c i ~ yj (2e//i)87ra 2 pt n ~ 600/s 
and ri ne i ~ 2Kn ~ 0.1/s, respectively. The favorable 
Tei/Finei ratio may enable evaporative cooling. A strin- 



gent constraint on the use of narrow- line optical Feshbach 
resonances is the proximity of the atomic transition. In 
the above example, the scattering rate of the PA laser 
photons due to the 1 So— 3 Pi atomic line is T s ~ 40/s 
~ r c i/15. Excessive one-atom photon scattering can 
cause heating and loss of atoms from the lattice. 

Another application of narrow line PA is potentially 
efficient production of cold molecules in the ground 
state. Our bound-bound Franck-Condon factor calcula- 
tions show that over 50% of the molecules electronically 
excited to the —222 MHz level decay to the last ground 
state vibrational level (distributed between J = 0, 2). 

In summary, we have measured a series of nine vibra- 
tional levels of the 88 Sr2 dimer near the 1 So— 3 Pi inter- 
combination transition. The 20 /is long lifetime of 3 Pi, 
combined with the magic wavelength optical lattice tech- 
nique, allowed direct probing of the least-bound state, 
and observation of thermal photoassociation line shapes 
even at ultracold /iK temperatures, with explicit dimen- 
sional effects. We have characterized the strengths of the 
molecular resonances, and shown that the least-bound 
state should allow a unique combination of extensive, yet 
low-loss optical control of the 88 Sr ground state scatter- 
ing length that may enable efficient evaporative cooling. 
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